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ABSTRACT 

We report results from deep observations ('^750 ks) of Tycho’s supernova remnant (SNR) with 
NuSTAR. Using these data, we produce narrow-band images over several energy bands to identify 
the regions producing the hardest X-rays and to search for radioactive decay line emission from “^^Ti. 

We find that the hardest (>10 keV) X-rays are concentrated in the southwest of Tycho, where recent 
Chandra observations have revealed high emissivity “stripes” associated with particles accelerated 
to the knee of the cosmic-ray spectrum. We do not find evidence of ^“^Ti, and we set limits on its 
presence and distribution within the SNR. These limits correspond to a upper-limit mass of 
M44 < 2.4 X 10~^Mq for a distance of 2.3 kpc. We perform spatially resolved spectroscopic analysis 
of sixty-six regions across Tycho. We map the best-fit rolloff frequency of the hard X-ray spectra, and 
we compare these results to measurements of the shock expansion and ambient density. We find that 
the highest energy electrons are accelerated at the lowest densities and in the fastest shocks, with 
a steep dependence of the roll-off frequency with shock velocity. Such a dependence is predicted by 
models where the maximum energy of accelerated electrons is limited by the age of the SNR rather 
than by synchrotron losses, but this scenario requires far lower magnetic field strengths than those 
derived from observations in Tycho. One way to reconcile these discrepant findings is through shock 
obliquity effects, and future observational work is necessary to explore the role of obliquity in the 
particle acceleration process. 

Subject headings: X-rays: ISM — ISM: supernova remnants — ISM: individual object: Tycho’s SNR 


1. INTRODUCTION 

Tycho’s supernova remnant (SNR G120.1-I-1.4; here¬ 
after Tycho) is widely believed to be the rem¬ 
nant of a Type la S N _ exp losion observ ed _ in 

1572 (l^tephenson fc GreenI 1200^ iBadenes et al.l 120061 : 
iRest et ^ I2OO8I : iKrause et al.l 120081) . Evidence of 
particle accelera t ion i n Tycho was found first by 
iPravdo fc SmithI ()1979H . who detected non-thermal X- 
ray emission up to 25 keV with HEAO-1. Initial Chan¬ 
dra X-ray Observatory images spatially resolved the non- 
thermal X-ray emission and showed it or iginates from 
narrow filaments around the rim of Tycho (iHwang et al.l 
I2OO2I : iBamba et al.ll2005t I Warren et al.ll2005D . Such fea- 
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tures are now thouKh t to be common among young SNRs 
(see [Reynold^ 1200^ for a review). A deep Chandra 
program revealed several non-thermal, high emissivity 
“stripes” in the projected interior of the SNR, a result 
which was interpreted as direct evidence of particles ac¬ 
celerated to the “knee” of the cosmic-ray ( CR) spec- 
tru m at around 3 Pe V (jEriksen et al.l I2011L although 
see iBvkov et al.l 120111 for an alternative interpretation). 
Recentl y. Tycho has been de tected in GeV y-rays with 
Fermi dGiordano et al.l 1201211 and in TeV y-rays with 
VERITAT iTAcciari et al.R^Oll . The y-ray spectrum is 


consistent with diffusive-shock acceler ation (DSA) either 
of GR proto ns (iBerezhk o et al. 1201311 or of two lepton 


populations l|Atovan fc Dermei]l2012ll . 


Recent X-ray studies of Tycho have reported the 
possible detection of radioactive d ecay lines of ^^Ti 
(|Troia et al.l 12014 IWang fe TI 1201 4 an d of Ti-K line 
emission at ^4.9 keV (|Miceli et al.ll2015[l . As the prop¬ 
erties of the '^^Ti (like yield, spatial distribution and 
velocities) in young SNRs probe directly the under- 
lying explosion mechan ism of the supernova (see e.g., 
iMagkotsios et all 1201011 , constraints on these parame¬ 
ters are useful to motivate and test sim ulations. The 
^^Ti isotope has a half-life of 60 years (jAhmad et al.l 
[Mil and originates from a-rich freeze-out near nu- 
clear statistical equilibrium (e.g., iTimmes et al.l 119961 : 
iWooslev et al.ll200l . Radioactive decay lines of ^'^Ti to 
■^"^Sc and then to are observable at ^68, 78, and 


core-collapse SNRs Cassiopeia A (e.g., 

Ivudin et al.il994 

IVink et ai.ll2001l:IRenaud et al.ll2006bl: 

Grefenstette et al.l 
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TABLE 1 

NuSTAR Observation Log 


ObsID 

Exposure 

UT Start Date 

40020001002 

339 ks 

2014 April 12 

40020011002 

147 ks 

2014 May 31 

40020001004 

262 ks 

2014 July 18 


120^ and SN 1987A (iGrebenev et al.ll2012HBoggs et alJ 
I2015r) . can also manifest itself through a 4.1 keV line 
which occurs from the filling of the inner-shell vacancy re¬ 
sulting from the decay of to "^'^Sc by electron capture. 
A feature at this energy has been detect ed in Chandra X- 
ray sp ectra of the young SNR Gl.9-1-0.3 (|Borkowski et al.l 
[20in[ [201111 . To date, no Type la SNRs yet have defini¬ 
tive detections o f the 68, 78, or 1157 keV lines associated 
with ^^Ti le.g.. rPupraz et al.l 11997 : llvudin et al.l 119991 : 

iRenaud et al.ll2006at iZoglauer et ani2015ll. 

In this paper, we present hard X-ray images and spa¬ 
tially resolved spectra from a set of NuSTAR observa¬ 
tions of Tycho totaling '^750 ks. Launched in 2012, NuS¬ 
TAR is the fi rst satellite to focus a t hard X-ray energies 
of 3-79 keV ([Harrison et al.ll2013f l. The primary scien¬ 
tific motivation of our observing program was to map 
and characterize the non-thermal X-ray emission in Ty¬ 
cho and to detect or constrain the spatial distribution of 
44Ti The paper is outlined as follows. In Section O we 
describ e the data reduction and analysis procedures. In 
Section [in we present the composite narrow-band NuS¬ 
TAR images of the hard X-ray emission in Tycho, and 
in Section we exploit NuSTAR images to set upper 
limits on the presence of '*^Ti. In Section (3.31 we report 
the results from a systematic spatially resolved spectro¬ 
scopic analysis across the SNR. In Section IH we discuss 
the implications of our results, specifically related to the 
particle acceleration properties of Tycho (in Section im 
and to "^^Ti searches in young SNRs (in Section iTlT) . We 
summarize our conclusions in Section [5l 

2. OBSERVATIONS AND DATA ANALYSIS 

Tycho was observed by NuSTAR three times from 
April-July 2014, as listed in Table [U with a total net 
integration of 748 ks. We reduced these data using 
the NuSTAR Data Analysis Software (NuSTARDAS) 
Version 1.3.1 and NuSTAR CALDB Version 20131223. 
We performed the standard pipeline data processing 
with nupipeline, with the stricter criteria for the pas¬ 
sages through the South Atlantic Anomaly (SAA) and 
the “tentacle”-like region near the SAA to reduce back¬ 
ground uncertainties. 

Using the resulting cleaned event files, we produced 
images of different energy bands using the ETOOL xse- 
lect and generated associated exposure maps using nuex- 
pomap. As Tycho is a bright, extended source, we opted 
to model the background and produce synthetic, energy- 
dependent background images for backgroun d subtrac¬ 
tion. We fo llowed the procedu r e deta iled in iWik et ahl 
([20141 1 and iGrefenstette et all ([20141 1 to estimate the 
background components and their spatial distribution. 
Subsequently, we combined the vignetting- and exposure- 
corrected FPMA and FPMB images from all epochs us¬ 
ing ximage. 

The combined images were then deconvolved by the 
on-axis NuSTAR point-spread function (PSF) using the 



Fig. 1. — Deconvolved, background-subtracted NuSTAR image 
from the combined <^750 ks observation of Tycho, where red is 
6—7 keV (which is dominated by Fe line emission) and blue is 10- 
20 keV. The SNR is >^9' in diameter. In this and all subsequent 
images, north is up, and east is left. 

maxJikelihood AstroLib IDL routine. The script employs 
Lucy-Richardson deconvolution, an iterative procedure 
to derive the maximum likelihood solution. We set the 
maximum number of iterations to 50, as more iterations 
did not lead to any significant changes in the resulting im¬ 
ages. We note that this routine assumes the data can be 
characterized by a Poisson distribution, but background 
subtraction causes the images to not follow strictly a 
Poisson distribution. Thus, the deconvolved images in 
Figures [T] and [5] are presented for qualitative purposes 
only, and we do not use them for any quantitative results. 

We performed spatially resolved spectroscopic analy¬ 
ses by extracting and modeling spectra from locations 
across the SNR. Using the nuproducts FTOOL, we ex¬ 
tracted source spectra and produced ancillary response 
files (ARFs) and redistribution matrix files (RMFs) from 
each of the three observations and both the A and B mod¬ 
ules (3 ObsIDS X 2 modules = 6 spectra per region). Fur- 
thermor e, we utilized the nuskybgd routines (presented in 
detail in iWik et al.l20l3) to simulate the associated back¬ 
ground spectra. We employed the addascaspec FTOOL 
to combine the six source pulse-height amplitude (PHA) 
files, ARFs, and background PHA files associated with 
each region. The six RM Fs were added using the ad drmf 
FTOOL. As discussed bv IGrefenstette et al. ([20I5D . the 
addition of spectra from different observations leads to 
systematic effects of order a few percent in the best-fit 
normalizations of the spectra. 

For comparison to the NuSTAR data, we analyzed 
archival Chandra and National Radio Astronomy Obser¬ 
vatory (NRAO) Very Large Array (VLA) observations 
of Tycho. The Chandra Advanced GGD Imaging Spec¬ 
trometer (AGIS) obse rved Tycho in April 2 009 for a total 
integration of 734 ks (lEriksen et aLllMTlIl . Reprocessed 
data were downloaded from the Chandra archive (ObsIDs 
I0093-I0097, I0902-I0906), and composite exposure- 
corrected images were produced of the 4-6 keV contin- 
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Fig. 2. — Deconvolved, background-subtracted NuSTAR images of three hard energy bands: 8-10 keV (left), 10—20 keV (middle), and 
20—40 keV (right). White contours on right panel are from the 10-20 keV band to guide the eye. 



Fig. 3. — Left: Three-color image, with 4—6 keV Chandra data in red (binned by a factor of eight), 8-10 keV deconvolved NuSTAR data 
in green, and 10—20 keV deconvolved NuSTAR data in blue. The hardest X-rays are brightest in the west of the S NR and coincide wit h 
the non-thermal X-ray “stripes” (denoted by the large red rectangle) seen in the deep Chandra data and analyzed bv lEriksen et aO II2011I1 . 
The small rectangle corresponds to the “hook” of non-thermal emission described in the text. Right: Three-color image, with 4-6 keV 
Chandra data in red, 10—20 keV deconvolved NuSTAR data in green, and 20-cm VLA data in blue. We note that the NuSTAR, Chandra, 
and VLA images have pixel sizes of 2.5", 3.9", and 4.2”, respectively. 


uum using the fluxjimage command in the Chandra In¬ 
teractive Analysis of Observations (ciAO) software Ver¬ 
sion 4.6. The VLA observed Tycho at 1.505 GHz on 10 
July 2001 in the C configuration for 174 minutes (pro¬ 
gram AR464, PI: Reynoso). The reduced image, with a 
beam size of 15.8" by 13.5", was downloaded from the 
NRAO/VLA image archive. 

3. RESULTS 
3.1. NuSTAR Imaging 

Figures [T] and [2] show the deconvolved, background- 
subtracted NuSTAR images of Tycho from the combined 
^^750 ks observation in several energy bands. NuS¬ 
TAR detects the “fluffy” Fe-rich eje cta found predom¬ 
inant l y in the northwest of Tycho (iDecourchelle et all 
120011 : (Vamaguchi et al.ll2014D as well as non-thermal X- 
rays with up to energies of ^40 keV. The hard (>8 keV) 
X-rays are predominantly distributed around the rim of 
the SNR, with the peak at the location of the “stripes” 

https://science.nrao.edu/facilities / vla/ar chi ve/index 


identified in Chandra observations (see Figure [3] for a 
comparison). 

In the right panel of Figure[3]we compare the 20-cm ra¬ 
dio morphology seen by the VLA to the 4-6 keV Chandra 
and 10-20 keV NuSTAR images. The radio emission is 
distributed in a ring interior to the narrow, non-thermal 
filaments at the SNR periphery found by Chandra. Gen¬ 
erally, the bright, hard X-ray knots detected by NuSTAR 
do not coincide with the radio emission. Overall, we find 
that the Chandra 4-6 keV band has a similar morphology 
as revealed in the NuSTAR images up to ^^20 keV. The 
hardest X-rays (20-40 keV) are concentrated around the 
stripes, as well as near a “hook” in the southwest of the 
SNR (denoted by the small red rectangle in Figure [3]). 

We also convolved the Chandra 4-6 keV image with a 
Gaussian to match the resolution of the NuSTAR data, 
and the comparison of the result to the 10-20 keV NuS¬ 
TAR data is shown in Figure SI As the 10-20 keV 
emission is almost exclusively non-thermal, the locations 
where the 4-6 keV emission is coincident reveals the lat¬ 
ter is substantially non-thermal in nature. The con- 
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Fig. 4. — Two-color image of Tycho to demonstrate the spatial 
coincidence of the 4-6 keV Chandra data (in red) and the 10-20 
keV deconvolved NuSTAR data (in green). We have convolved the 
Chandra image with a Gaussian function to match the resolution 
of the NuSTAR image. As the emission in the 10—20 keV band 
is almost exclusively non-thermal, locations where the 4-6 keV 
emission is coincident indicates that the latter is also non-thermal 
in nature. 

verse scenario (where 4-6 keV emission is found with- 
out 10-20 keV emission) indicates the softer band may 
be dominated by thermal emission. We find that the 
4-6 keV emission is largely coincident with the harder 
band, except in the northwest, where previous X-ray 
studies have found significant emission from th e ejecta 
(|Decourchelle et aT1l200HlYamaguchi et al.l[2014ll . Thus, 
we conclude that the 4-6 keV Chandra band is likely 
dominated by non-thermal across Tycho except in the 
northwest. 

In Figure [Sj we compare the NuSTAR 20-40 keV 
emission to the GeV and TeV centroids rep orted by 
IGiordano et al.l (|2012ll and lAcciari et al.l (|2011ll . respec¬ 
tively. The GeV localization (from Fermi observations 
with a 95% confidence level) coincides with the hard X- 
rays from the “hook” feature in the southeast of the SNR, 
although we note that the “stripe” region is just outside 
the 95% error circle. The TeV detection is consistent 
with a point source at the 0.11° 6.6' PSF of VERI¬ 

TAS, which is of the same order as the angular extent 
(^8') of the SNR. Thus, it is not possible with the cur¬ 
rently reported VERITAS observations to determine the 
origin of the TeV emission. 

3.2. Upper Limits on Ti 

We do not find any point-like or extended emission in 
the remnant using the 65-70 keV band containing the 68 
keV '^^Ti emission line. We se arched using the method 
from iGrefenstette et all (|20I4ll , where we compared the 
observed image to a simulated background image. We 
searched on multiple spatial scales by convolving both 
the NuSTAR image and the simulated image with a top 
hat function, varying the radius of the function from 3-15 
NuSTAR pixels (roughly 7.5 to 30 arcseconds). We com¬ 
puted the Poisson probability that the observed counts 
in each pixel could have been produced by the predicted 
background flux, where a low probability indicates that 



Fig. 5.— Comparison of the high-energy detections of Tycho, 
with the Fermi 95% confidence error circle denoted in red (for 
photons of energy E > 1 GeV; from IGiordano et al.l 120121 ) , the 
VERITAS error ci rcle in green (for energies E > 1 TeV; from 
lAcciari et ahllTOllI) . and the NuSTAR 20—40 keV band contours in 
blue. The image is the background-subtracted NuSTAR 20—40 keV 
image smoothed with a Gaussian function of width cr = 5 pixels. 

the pixel likely contains a source count. We do not find 
any evidence for excess emission on any spatial scales. 

We place upper limits on the diffuse extended emis¬ 
sion in Tycho arising from the 68 keV line. We extracted 
spectra and produced an ARE and RMF for source re¬ 
gions with radii from I to 5 arcminutes centered on the 
remnant. For example, in Figure [6l we give the spectra 
and background model of the FPMA and FPMB data 
for a 4' radius extraction region. For each radius, we 
modeled the source using a combination of a power-law 
continuum with two Gaussian lines of fixed relative in¬ 
tensity centered at 68 and 78 keV and fit the data over 
the 10-79 keV bandpass. We do not find any evidence 
for line emission in any of the source regions. 

Our upper limits are dependent on the intrinsic width 
of the "‘^Ti lines. We computed the 3-tT confidence limits 
on the line normalizations for fixed line widths of O.I, 
I, or 3 keV, correspond to velocities of «733, 7330, and 
22,000 km s“^, respectively. The resulting 3-a upper 
limits as a function of radius are shown in Figure [71 In 
all cases, the lower limits are zero. 

To effectively compare these NuSTAR upper limits 
with results from other missions, we convert all previ¬ 
ous measurements and limits to epoch 2014. In Ta¬ 
ble m we list these values for observations with the 
Compton Gamma-ray Observatory/COMPTEL, INTE- 
GRAL/IBIS, and Swift/JiAT. To obtain epoch 2014 
fluxes (which are a lso given in Table [2]) , we adopt a half- 
life of 58.9 years ([Ahmad et al.l [200a) . The NuSTAR 
upper limits may be used in conjunction with the previ¬ 
ous results to assess the likely spatial distribution of the 
■^^Ti. In particular, our upper limits for radii ^ 2' are less 
than the reported Swift/BAT detection (at epoch 2014) 
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Fig. 6.— Top Left. Integrated FPMA source spectra (black) and background model (red) of a 4' radius region. Top right. FPMA 
background-subtracted spectra. Right panels: Same as left panels for FPMB. Vertical red dashed line shows the rest-frame energy of the 
67.86 keV Ti line. We do not find any excess emission at this line energy relative to the background model, and thus we set upper limits 
on 44Ti. 


of (1.2±0.6)xl0“^ ph cm“^ s“^, suggesting that the ‘^'^Ti 
is distribut ed over a larger rad ius of Tycho. 

Recently. iMiceli et al.l (|2015D reported detection of the 
Ti-K line complex at ~4.9 keV using stacked XMM- 
Newton spectra from Tycho. In particular, these authors 
identified the Ti feature in two regions of enhanced emis¬ 
sion from Fe-group elements, in the northwest and in the 
southeast of the S NR (regions #1 and #3 in Figure 2 of 
IMiceli et al. 1120151) . Motivated by this result, we searched 
for the 68 keV line at these locations in the NuSTAR data 
(the region files were provided graciously by M. Miceli). 
We did not detect any signal above the background from 
either region, with lower limits of zero in all cases. The 3- 
tr confidence flux upper limits are 3.9x10“® ph cm“^ s“^ 
and 7.3x10“^ ph cm“^ s“^ for regions #1 and #3, re¬ 
spectively, assuming a line-width of 1 keV. 

3.3. Spatially Resolved Spectroscopy with NuSTAR 

We performed spatially resolved spectroscopic analyses 
by extracting and modeling spectra from sixty-six 1' by 
1' boxes across the SNR (see Figured]). To evaluate the 
non-thermal emission, we fit the spectra in the 10-50 
keV range with an absorbed srcut model. In our fits, we 
have fixed the absor bing column to iVn = 7 x 10^^ cm“^ 
(|Eriksen et al.ll20Tli) . and we have grouped the spectra 



Fig. 7.— The black, red, and green points show the 3-cr confi¬ 
dence limits as a function of radius for line widths a of 0.1, 1, and 
3 keV (corresponding to velocities of ~730, 7300, and 22,000 km 
s-l). 

to a minimum of 30 counts bin“^. 

The srcut model describes the spectrum as that 
radiated by a power-law energy distribution of elec- 
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TABLE 2 

Previous Measurements and Limits'^ on 


Instruments 

Years of 
Obs. 

Mean Year 
of Obs. 

Flux at Mean Year 
(xl0“^ ph cm“^ s“^) 

Flux at Epoch 2014 
(xl0“® ph cm~^ s“l) 

CGRO/COMPTEL 

1991-1997 

1994 

<1.6 

<1.3 

INTEGRAL/IBIS 

2002-2006 

2004 

<1.5 

<1.4 

Swift/BAT 

2004-2013 

2008.5 

1.3±0.6 

1.2±0.6 


® CGiJO/COMPTEL upper limit is 2-cr, from|Dupraz_et_ay II1997I1 a ndllvudin et alJ II1999I) : INTE¬ 
GRAL/IBIS upper limit is S-cr, from lHe naud et al.l Il20 06all ari d lWang Hz. Lil II2014I1 : the Swift/BAT 
detection is the 90% confidence limit. from lTroiTetaTT (1201411 . 



Fig. 8.— Background-subtracted, merged 3-20 keV image of Ty¬ 
cho with 1' by 1' boxes overplotted denoting the regions where 
NuSTAR spectra were extracted and modeled. 

trons with an exponentia l cutoff at energy Emax 
(|Revnolds fc Keohand Il999ri . That spectrum cuts off 
roughly as exp {— \/v/VmWos) (though the fitting uses 
the complete expression). The rolloff photon energy 
Eroiioff = ^t^roiioff IS related to the maximum energy of 
the accelerated electrons Emax by 


E 


max 


120 


^t^rollofF 

1 keV 



TeV 


( 1 ) 


where B is the magnetic field strength and we 
have corrected for a s mall numerical error in 
[Reynolds fc Keohand (|1999D . 

In XSPEC, the srcut model is characterized by three 
parameters: the rolloff frequency iZroiioff, the mean radio- 
to-X-ray spectral index a, and the 1 GHz radio flux den¬ 
sity Eighz- To limit the free parameters in the fit, we 
estimated Eiqhz in each of the sixty-six regions by mea¬ 
suring the 1.505 GHz flux density in the NRAO/VLA 
data assuming a r adio spectral index of a = —0.65 
(iKothes et al.ll200^ . The derived flux densities at 1 GHz 
using this procedure are listed in Table |31 We note that 
without additional single-dish observations, the interfero¬ 
metric VLA data are missing flux from the largest scales. 
Based on previous radio studies of Tycho, we estimate 
the missing flux is of order ~10% of the values listed in 
Tabled 

Table [3] lists the spectroscopic results, giving the best- 


fit rolloff frequency i^roiioff and the and degrees of 
freedom (d.o.f.) for all sixty-six regions. We find that 
all the regions are fit well by the srcut models, with 
X^/d.o.f. = ~0.82-1.26 with ^100-200 degrees of 

freedom per region. Eigure [9] shows example spectra 
from six of the regions, and Eigure 13.31 maps the best- 
fit rolloff frequency iZroiioff values across Tycho. We find 
substantial variation of a factor of five in lymiioS- The 
largest values of t'roUoS ~ 3 x 10^’’ Hz (corresponding to 
a rolloff energy Aroiioff =1.2 4 keV) are concentra ted in 
the west of the SNR where lEriksen et al.l (|201ir ) iden¬ 
tified the “stripes” in Tycho. Furthermore, we find el¬ 
evated rolloff frequencies in the northeast of the SNR, 
with t'roUoff ~ 2 X 10^^ Hz (corresponding to E^oWoS ~ 
0.83 keV), where the SN R is thought to be interacting 
with molecular material (iR.evnoso et ahllTo^ iLee et al.l 
1200411 . The lowest values of t'roiioff 10^^ Hz (corre¬ 
sponding to Eroiioff ~ 0.42 keV) are found in the south¬ 
east and north of the SNR. For the range of i?roUoff ~0.4- 
2.0 keV we find in our spatially resolved spectroscopic 
analyses. Equation [T] gives a maximum energy Amax ~ 
(75 — 170)(H/^G)“^/^ Te V. Assuming a mag netic field 
B = 200/rG across Tycho (|Parizot et al.ll200W . this rela¬ 
tion gives Emax = 5-12 TeV. 

Previous X-ray studies of Tycho have also mea¬ 
sured t'roiioff over the entire SNR or along prominent 
sync hrotron filaments. Using ASCA data from Ty¬ 
cho, iRevnolds fc Keohand (|1999fl found t'roiioff ~ 8.8 x 
10^® Hz, and using the h rst Chandra observation of Ty¬ 
cho, iHwang et aLl (|2002ll derived even lower t'roiioff to¬ 
ward the northeast, northwest, and southwest rim, with 
^'rolloff ~ 3-7 xlO^® Hz. More recently, in t heir analysis of 
the de ep Chandra observations of Tycho, lEriksen et al.l 
(|2011f) fit absorbed srcut models to the hard “stripes” 
and a series of filaments projected in the south of Tycho 
(corresponding to our regions ^51 and 52) which they 
called the “faint ensemble”. They found the stripes have 
t'roiioff 2 X 10^® Hz while the faint ensemble had softer 
spectra, with t'roiioff ~ 3 x 10^^ Hz. 

Our NuSTAR results confirm the stripes have harder 
emission than other regions in Tycho, but our 90% con¬ 
fidence limits exclude t'roiioff ~ 5 x 10^^ Hz there. The 
softest regions (with 10^^ Hz) are co i npara - 

ble to the values found bv iR.evnolds fc Keohand (Il999f l 
integrated over Tycho. The discrepancy between our 
90% confidence limi ts and the hard stripes analyzed by 
lEriksen et al.l (l2011h may be due to two issues. First, 
the angular resolution of NuSTAR is not sufficient to re¬ 
solve the individual stripes, and thus our spectra appear 
softer as they are averaging ov er larger areas of the SNR. 
Secondly. lEriksen et al.l (j2011h fit the Chandra data from 
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TABLE 3 

Spatially Resolved Spectroscopy Fit Results Using SRCUT^ 


Region 

^IGHz 

(Jy) 

Urolloff 

(xlO^'^ Hz) 

X^/dof 

Region 

i^lGHz 

(Jy) 

Urolloff 

(xlOi'^ Hz) 

X^/d.o.f. 

1 

0.55 

1.91±0.06 

109/126 

34 

0.42 

1.55±0.05 

117/110 

2 

1.32 

1.17±0.03 

166/138 

35 

0.53 

1.59±0.05 

138/123 

3 

1.29 

1.08±0.03 

130/122 

36 

0.75 

1.63±0.04 

133/147 

4 

0.89 

1.22±0.04 

116/117 

37 

1.18 

2.25±0.04 

247/213 

5 

0.97 

1.33±0.04 

129/128 

38 

0.58 

3.43±0.08 

204/186 

6 

1.85 

1.17±0.02 

156/158 

39 

1.38 

1.11±0.03 

161/132 

7 

1.51 

1.21±0.02 

161/154 

40 

0.90 

1.20±0.03 

135/121 

8 

0.97 

1.40±0.03 

163/140 

41 

0.74 

1.66±0.04 

164/142 

9 

1.27 

1.30±0.03 

152/148 

42 

0.76 

1.41±0.03 

133/127 

10 

1.39 

1.33±0.03 

144/143 

43 

0.46 

1.62±0.05 

137/114 

11 

0.40 

2.01±0.08 

93/108 

44 

0.48 

1.78±0.05 

132/121 

12 

0.34 

2.11±0.08 

134/109 

45 

0.85 

1.78±0.04 

155/159 

13 

1.51 

1.09±0.02 

145/141 

46 

1.08 

2.59±0.04 

208/218 

14 

1.05 

1.29±0.03 

144/135 

47 

0.47 

3.66±0.09 

191/172 

15 

0.74 

1.45±0.04 

130/128 

48 

0.87 

1.22±0.04 

133/117 

16 

0.55 

1.65±0.05 

110/120 

49 

1.05 

1.14±0.03 

134/123 

17 

0.63 

1.83±0.05 

124/134 

50 

1.16 

1.36±0.03 

154/154 

18 

1.15 

1.48±0.03 

162/152 

51 

1.03 

1.34±0.03 

154/147 

19 

1.20 

1.56±0.03 

152/159 

52 

0.69 

1.58±0.04 

147/135 

20 

0.14 

5.38±0.27 

145/115 

53 

0.71 

1.81±0.04 

140/146 

21 

0.67 

2.52±0.06 

147/160 

54 

1.01 

2.25±0.04 

217/198 

22 

1.38 

1.13±0.02 

165/136 

55 

1.32 

2.16±0.03 

214/205 

23 

0.85 

1.25±0.03 

136/120 

56 

0.19 

4.22±0.18 

110/121 

24 

0.62 

1.42±0.04 

117/116 

57 

1.15 

1.33±0.03 

163/144 

25 

0.54 

1.47±0.04 

105/114 

58 

1.40 

1.27±0.02 

164/158 

26 

0.58 

1.58±0.04 

122/125 

59 

1.36 

1.37±0.03 

157/153 

27 

0.80 

1.54±0.04 

120/143 

60 

1.08 

1.57±0.03 

157/157 

28 

1.21 

1.58±0.03 

162/167 

61 

1.04 

1.74±0.04 

127/155 

29 

0.54 

2.64±0.07 

179/153 

62 

0.48 

2.45±0.07 

149/139 

30 

1.09 

1.55±0.03 

188/151 

63 

0.49 

1.71±0.06 

130/115 

31 

0.97 

1.18±0.03 

131/118 

64 

0.76 

1.60±0.04 

126/133 

32 

0.66 

1.39±0.04 

113/117 

65 

0.53 

1.69±0.05 

127/117 

33 

0.64 

1.37±0.04 

135/116 

66 

0.09 

4.63±0.31 

101/94 


^ Error bars on the rolloff frequency Ufoiioff represent the 90% confidence range. 







Fig. 9. — Example spectra and fits using an absorbed srcut model for six regions. Although the spectra were fit over the range 10-50 keV 
where the non-thermal emission dominates, the spectra above span 3—50 keV to demonstrate the adequacy of the srcut model at predicting 
the non-thermal continuum across this energy range. 
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Fig. 10.— Map of roll-off frequency i^rolloff across Tycho, with 
black contours from the 3-20 keV image in Figure [S] The softest 
spectra (i.e., those with the smallest best-fit values of i.'rollofi) Eirs 
found in the southeast and north of Tycho, while the hardest spec¬ 
tra are in the northeast and the west. The hardest spectra (with 
^^rolloff ~ 3 X 10*^^ Hz) are coincident with the “stripes” found in 
deep Chandra observations of Tycho bv lEriksen et aO II2011I) . 


4.2-10 keV, and these spectra may not have had reliable 
leverage to model the high-energy rolloff. Additionally, 
we note that the effective area of Chandra drops precip¬ 
itously above 8 keV. 

We note that we have assumed a constant radio 
spectral index of a = —0.65, the value derived by 
IKothes et al.l (120061) based on the integrated flux den¬ 
sities at 408 MHz and 1420 MHz. However, the radio 
spectrum is actuall y concave, and thus, a ma y flatten at 
higher frequencies (iRevnolds fc Ellisoi] 119921) . If we in¬ 
stead adopted a = —0.5 in our analysis, the best-fit val¬ 
ues of i^roiioff change substantially, decreasing by about a 
factor of five. However, the assumption of a = —0.5 is 
not favored statistically based on the computed from 
those fits. For example, region ^21 has = 235 (147) 
with 160 degrees of freedom in the a = —0.5 (—0.65) 
case. Additionally, the a = —0.5 fits dramatically over¬ 
predict the continuum flux below 10 keV, whereas the 
a = —0.65 fits adequately describe the continuum down 
to 3 keV (see Figure 0 ). Generally, regardless of the 
choice of a, the qualitative trend in the spatial distribu¬ 
tion of r'roiioff (e.g., with greater r'roiioff in the west of the 
SNR) remains the same. 

4. DISCUSSION 

4.1. The Detection of'^^Ti in Young SNRs 

The yield, spatial distribution, and velocity distribu¬ 
tion of ^^Ti in a young SNR are a direct probe of the 
underlying explos ion mechanism of the originating su¬ 
pernova (see e.g.. iMagkotsios et al.ll201(3l) . Among core¬ 
collapse SNRs, Cassiopeia A remains the only Galactic 
source with a conhrmed detect ion of ^^Ti (jivudin et al.l 
11994 iGrefenstette et all I2014D . and SN 1987A in the 
Larg e Magellanic Cloud a l so has a robust ^ ^Ti detec¬ 
tion (|Grebenev et al.l[2012l : iBoggs et al.ll2015D . To date. 


the only probabl e detection of ^"^Ti in a Type la SNR 
was reported by iTroia etldl (|2014D . who found a 4 -(t 
level excess above the continuum in t he 60-85 keV band 
in Swiftl^KY observations of Tycho. iWang fcTil (j2ni4[ ) 
also found a marginal 2.6 -(t detection of excess emission 
in the 60-90 keV band using INTEGRAL/lYilS^ although 
their 3 -(t upper limit is con sistent with the prev ious IN¬ 
TEGRAL studies of Tycho (jRenaud et al.l[200^ . Other 
searches for "^“^Ti in Type la SN Rs have yielded upper 
limits, such as the recent work bv IZoglauer et al.l (I2015D 
using NuSTAR observations of the young SNR Gl.9-1-0.3. 

Based on our NuSTAR analyses, we do not detect '^‘^Ti 
in Tycho. We looked both for s mall scale “clumpy” emis- 
sion as found in Cassiopeia A (IGrefenstette et ani2014D 
as well as for diffuse emission throughout the SNR. We 
Hnd upper limits which are more st r ict than the previ- 
ous COMPTEL (iDupraz et al.lll99^ llyudin et al.l 1999 ) 
and /ATEGRAL (|Renaud et al.ll2006al : IWang fc Lil2014 ) 
values, particularly if the ^^Ti is concentrated in the 
central 2 ' and/or has slow or moderate velocities (i.e., 
the 0.1 keV and 1.0 keV line width cases). For exam¬ 
ple, our 3-cr limit in the 0.1 keV (1.0 keV) case at an 
enclosed radius of 2' is 5.0x10“® (1.0x10“®) ph cm“^ 
s“^. To be consis tent with the poss ible detection of ^"^Ti 
with Swift/BAT (|Troia et al.ll20l4) . our results indicate 
that the ^^Ti is expanding at moderate-to-high velocities 
(7000 km s“^) and/or is distributed across the majority 
(>3') of the SNR. 

We have also fo und no ^^Ti in the regions where 
iMiceli et al.l (I2015D reported detection of the Ti-K line 
complex using XMM-Newton. This ~4.9 keV feature 
arises from shock-heated, stable "^®Ti and ®®Ti. Thus, the 
lack of a co-located detection of ^"'^Ti may indicate that 
the radioactive isotope does not trace the Fe-rich ejecta. 
Alternatively, the divergent results may suggest that the 
4.9 keV feature arises from transitions of Ca xix and 
Ca XX rather than from stable Ti. In the future, these 
interpretations can be tested via deeper NuSTAR ob¬ 
servations to reveal the spatial distribution of '*^Ti, and 
through high-resolution X-ray spectroscopy (with e.g., 
Astro-H) which may distinguish between the Ga and Ti 
spectral features. 

The upper limit flux we have measured of the 67.9 keV 
line can be employed to set an upper limit on the "^^Ti 
yield in the SN explosion. The '^‘^Ti mass M 44 is re- 
lated to the flux Fe, H in the 67.9 keV ^"^Ti line (e.g., 
iGrebenev et al.|[20l3) by 

M44 = 47 rfi^( 44 mp) r exp(t/r) Fes ( 2 ) 

where d is the distance to the source, nip is the mass 
of a proton, r is the mean lifetime of the radioac¬ 
tive titanium (which is related to the half-life ti /2 by 
r = ti/ 2 /ln 2 w 85 years for ^“^Ti: I Ahmad et al.ll2(I(3^ . t 
is the age of the source, and W is the emission efficiency 
(average numbers of photons per decay) for a given emis¬ 
sion line (W = 0.93 for the 67.9 keV ^"^Ti line). Assuming 
a distance D = 2.3 kpc and an age oi t = 442 years for 
Tycho, an upper limit of Fes < 2 x 10“® ph cm“^ s“^ 
corresponds to M44 < 2.4 x 1 O“"‘M 0 . Here, we adopted 
the upper limit Fqs in the case where the ^“^Ti is located 
within a radius of 3' and has a line width of 1 keV, as 
these most closely match the expected scale and veloc- 
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ity of the shell. For the range of distance estimates in 
the literature oi D = 1.7-5.0 kpc to Tycho, we find up¬ 
per limits of M 44 < 1.3 X and M 44 < 10~^Mq, 

respectively. 

These upper limits can be compared to the predictions 
of ^"^Ti synthesized in Type la SN models. Generally, 
the amount of "‘'^Ti produced in a Type la SN depends 
on both the progenitor systems as well as the ignition 
processes. Simulations predict ^^Ti masses ranging from 
in centrally ignited pure deflagration models 
and ^ 6 X 10~^Mq for off-center delayed detonation mod¬ 
els (e.g., Ilwamoto et al.l iMaeda et al.l l2010l) . In 

sub-Chandrasekhar models, greater "‘'^Ti yie Ms are ex- 
pected, with estimates of > 1 O“^M 0 (e.g., iFink et ^ 
l201Ct iWooslev fc KaseiJ 120111: iMoll fc WooslevI 1201311 . 
The mass limits estimated above are most consistent with 
the models that produce small or moderate yields of ^'^Ti 
and disfavor the sub-Chandrasekhar models. 

4.2. Particle Acceleration in Tycho 

To explore the relationship between the non-thermal 
X-ray emission detected by NuSTAR and shock proper¬ 
ties, we compa re th e spatially resolved spectroscopic re¬ 
sults of Section 1531 to measurements of the shock expan¬ 
sion and the ambient density rio in Tycho from previous 
studies. 

Expansion measurements of Tycho’s rim have been 
done at radio, optical, and X-ray wavelengths. Us- 
ing 1375-MHz VLA data over a 10-year baseline, 
iRevnoso et al.l (|1997ll measured expansion rates of ~0.1- 
0.5" yr“^, depending on the azimuthal angle. These 
authors found the slowest expansion occurs in the 
east/northeast, wh ere the forward shock is interactin g 
with a dense cloud (|Revnoso et aLlHo^lLee et al.ll2004ll . 
and the fastest expansion is along the western rim. Quan- 
titatively similar results w e re obt ained in the optical by 
Kamper fc van den BerghI (1197811 and in the X-ray by 
Katsuda et al.l (|2oTo[ l. Using the measured expansion in - 
dex and SNR evolutionary models. iKatsuda et al.l (1201011 
inferred a pre-shock ambient density of Uolz 0.2 cm“^, 
with probable local variations of factors of a few based on 
the substantial differences in the proper motion around 
the rim. 

The shock expansion is described by the radius of the 
shock with time, i? oc t™, where m is the expansion pa¬ 
rameter. From this expression, the shock velocity is re¬ 
lated to m by Us = mR/t. Thus, the m is a useful indica¬ 
tor of shock velocity, yet it does not require assumptions 
about the distance to the source. Generally, m is close 
to unity when the expansion is undecelerated (as in e.g., 
the free expansion phase of a young SNR). As the sho ck 
undergoes deceleration, m ^0.6-0 .8 (|ChevalieHll982al lbll. 
and m = 0.4 once the shock reaches the Sedov phase. 
Finally, as the shock becomes radiative, m drops below 
0.4. While interior pressure is still important, m = 0.3 3 
(“pressure-driven snowplow”; e.g. iBlondin et aLllIO^ . 
and then falls to m = 0.25 when the s hock transitions t o 
a pure momentum-c onserving stage (iCiofR et al.l [19^ . 
In the case of Tvcho. IRevnoso et al.l (119971) demonstrated 
the expansion parameter m varies around the rim from 
m ~ 0.25 in the east (where the shock is decelerating 
due to the collision with a dense clump known as knot q : 
[Kamper fc van den Berghlll978l : iGhavamian et al.ll200C)i) 



Fig. 11.— Rolloff energy F'rolloff (determined in Section 13.311 
of regions around the rim of Tycho versus exp ansion parameter 
m (as listed in Table 2 of IRevnoso et aT1ll997l) . Overplotted are 
three curves which represent different dependences of F^rolloff on m\ 
-^rolloff (cyan solid line), F'rolloff (red dashed line), and 

-^rolloff (blue dotted line). These dependences are derived 

using the expressions in Equations fblfS] and given that R cc , so 
Vg = mRIt oc m. The former two trends are those expected if the 
maximum energy of the electrons is limited by the age of the SNR 
or by the timescale of the radiative losses, respectively. The latter 
curve represents the loss-limited case if the B-field strength is not 
constant (and instead B oc Ug oc m). The Eroiioff trend does 

not match the data for large m, and thus, the age-limited scenario 
is more consistent with our results. 

up to TO 0.8 in the southwest. 

In Figure [Til we show the rolloff energy (found in Sec- 
tion 13.311 versu s the e xpansion parameter m obtained by 
IRevnoso et al.1 (|1997[) using VLA data. To produce this 
plot, we have considered only our regions around the 
rim of Tycho, and we have adopted the mean of the 
expansion parameters over the azimuthal angles which 
correspond to those regions using the data in Table 2 of 
IRevnoso et al.lll9^ At to > 0.4, the rolloff energy ap¬ 
pears to increase with to, suggesting the least decelerated 
(fastest) shocks of Tycho are the most efficient at accel¬ 
erating particles to high energies. The minimum rolloff 
energy is around to = 0.4 where the SNR has reached its 
Sedov phase. At to = 0.25 where the shock is encounter¬ 
ing dense material, the SNR is slightly more efficient at 
accelerating particles, reaching half of the rolloff energy 
of the high-TO regions. 

In Figure da we plot the rolloff energy Eroiioff ver¬ 
sus the post-shock density n^, (left panel) and versus the 
shock velocity Vs (r ight panel). B o th of these quanti¬ 
ties were derived bv lWilliams et al.l (I2013D using Spitzer 
Space Telescope 24 and 70 fim d ata as well as t he ra - 
dio and X-ray proper mo tions of IRevnoso et al.l ()1997[ 1 
and IKatsuda et al.l ()2010[ ) . We find that the regions of 
large Aroiioff have the lowest post-shock densities (with 
77,0^0.25 cm“^) and the largest shock velocities (with 
Vs it 3500 km s“^). The two outliers to these trends 
(those with uq >1 cm“^, Vs < 2200 km s“^, and 
Aroiioff >0.8 keV) correspond to regions #12 and 21, 
where the shock is interacting with knot g. 
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Fig. 12.— Rolloff energy ifrolloff versus post-shock density rio (left) and shock velocity Vs (right) derived bv l Williams et II2013I 1 adopting 
a distance of D = 2.3 kpc. The velocities scale linearly with the assumed distance, but the trend of E^aWoS with Vs is independent of 
distance. Overplotted on the right panel are three curves which represent different dependences of ifrollofl on shock velocity Vg: ifrolloff 
(cyan solid line), ifrollofl oc Vg (red dashed line), and Frolloff nj (blue dotted line). The former two trends are those expected if the 
maximum energy of the electrons is limited by the age of the SNR or by the timescale of the radiative losses, respectively. The latter curve 
represents the loss-limited case if the B-field strength is not constant (and instead B oc Us). The Frolloff trend does not match the 

data, and the age-limited scenario is more consistent with our results. We note that the three no >1 cm“® points correspond to regions 
#12, 10, and 21 (in order of increasing no). The two Ug 5 2200 km s“^ points correspond also to regions #12 and 21, which are located 
where the forward shock is interacting with the dense clump know as knot g IKamper fc van den Berghlll979 : IGhavamian et J]|2000ll . 


We note that the conversion of proper motions to shock 
velocities requires an assumption about the distance to 


erature ranging from 1.7- 

-5.0 kpc (e.g.. lAlbinson et al. 

1986 

:ISchwarz et al.lll995l 

Volk et al.ll2008HHavato et al. 

201C 

; iTian & LeahvI 12011 

; iSlane et al.l 2014D. For the 


ployed a distance of D = 2.3 kpc, as suggested by 
IChevalier et al.l (I1980D . The shock velocities scale lin¬ 
early with the assumed distance, so the choice of dis¬ 
tance will shift the points in the right panel of Figure [T^] 
accordingly. 

The relationship between Eroiioff and Vs in Figure [T^] 
can be used to examine the mechanism limiting the maxi¬ 
mum energy of electrons Emax undergoing diffusive shock 
acceleration. The timescale tacc of the acceleration is 
set by the lesser of three relevant processes: the age of 
the SNR (i.e., the finite time the electrons have had to 
accelerate), the timescale of the radiative losses of the 
electrons, and the timescale for the particles to escape 
upstream of the shock. In each case, Emax (and con¬ 
sequently, Eroiioff) h as different dep endences on physi¬ 
cal parameters (from iRevnoldsI 120081 Equations 26-28), 
specifically Vs'. 


Fniax(age) oc VgtB{qRj) ^ 

(3) 

£'i„ax(loss) oc VsiqRjB)-^/"^ 

(4) 

F'niax(6Sc) OC 7?Amax- 

(5) 


In the above relations, t is the age of the SNR, -q is the 


gyrofactor (where q = 1 is Rohm diffusion and ?7 >1 is 
“Bohm-like” diffusion), and Amax is the maximum wave¬ 
length of MHD waves present to scatter electrons. Rj 
parametrizes the obliquity dependence of the accelera¬ 
tion: given that the shock obliquity angle 0Bn is the 
angle between the mean upstream B-field direction and 
the shock velocity, then the quantity Rj is defined as 

Rj — ^acc(^Bn)/^acc(^Bn — 0)- SiuCe Arolloff CX! j 

the above relations can be rewritten in terms of Eroiios: 


Eroiioeiage) (X Vgt'^B^{qRj) ^ 

(6) 

Froiioff(loss) oc vl{qRj)~^ 

(7) 

AfollofF(asc) oc B Amax- 

(8) 


In the right panel of Figure [121 we plot curves for 
-E'roiioff oc Vg (solid light blue line) and Epoiioff oc 
(dashed red line), for comparison. Additionally, since 
Vs = mR/t oc m (where R is the radius of the shock and 
t is the age of the SNR), we have plotted the same curves 
on Figure [TT] as well. The two dependencies of Froiioff on 
Vs reflect those of the loss-limited and age-limited rela¬ 
tions above, assuming constant B and Rj. We find that 
the oc Wg curve is too shallow to match the high-velocity 
{vs ^ 3500 km s“^) points, and the oc curve is a better 
descriptor of the data. We note that if B is not constant 
and is instead amplified by a constant fraction of the 
shock kinetic energy, then B ocvs, and Equation [5] gives 
Aroiioff cc vl- We have plotted this steeper curve in Fig¬ 
ures dn and [12] as well, and it adequately describes the 
data. Due to the relatively large error bars on the shock 
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velocity Vs, we are not able to distinguish which curve 
(Aroiioff oc 'I’s Aroiioff oc vj) is better statistically, but 
both are much more successful than the loss-limited case 
(with Aroiioff oc fg: with no dependence on magnetic-field 
strength B) and the escape-limited case (with no Aroiioff 
dependence on Vs). Thus, our spatially resolved spectro¬ 
scopic results are most consistent with the age-limited 
scenario, whereby the maximum electron energy Amax is 
determined by the finite time the electrons have had to 
accelerate. In this case, the accelerated electrons achieve 
the same maximum energy as the protons/ions because 
the electrons are not limited by their radiative losses. 

However, previous studies of Tycho have argued the 
SNR is lik ely loss-limited instead of age-limited (al- 
though see [Reynolds fc Keohanel 119991) . In particular, 
IParizot et al.l (|2006[1 and iHelder et al.l (|2012li calculate 
the synchrotron loss timescale and found it to be a few 
percent of the age of Tycho. Furthermor e , rece nt hy- 
drodynamical simulations by iSlane et al.l (l2014f l show 
the loss-limited scenario is consistent with the observed 
broadband spectrum of Tycho. We note that these mod¬ 
els were matched to the global spectra across Tycho as 
the gamma-ray emission is not resolved. Thus, local vari¬ 
ations in the gamma-ray emission is possible, with the 
implication that the electrons and protons may achieve 
the same maximum energies in certain locations. 

For the acceleration to be age limited, the synchrotron 
loss time tsyn for X-ray emitting electrons must be larger 
than the remnant age, t =442 years for Tycho. As a con¬ 
sistency check, we estimate the magnetic field strength 
B necessary for tsyn > t and compare it to that measured 
observationally. The synchrotron loss time Avn is relat ed 
to the Aniax and Aroiioff by (adapted from I Vin 3 [Ml) 


^syn — 12.5 


tsyn = 104 


E„ 


100 TeV 


)■ 


B 


Arolloff \ 

1 keV J 


- 1/2 


100 nG 
B 

100 /iG 


-2 


yr 


(9) 


-3/2 


yr (10) 


Adopting the highest rolloff frequency from our spectral 
analysis, t'roiioff = 5.38 x 10^^ Hz (from region #20), we 
derive B = 29 ijlG. This value is far below the mea¬ 
surements bv IParizot et al.l (l2006l) oi B = 200/rG based 
on the assumption that the X-ray filament widths are 
limited by the synchrotron loss time. However, it is pos¬ 
sible that filament width s are governed instead by mag¬ 
netic field damping (e.g. , iPohl et iil]l2005l : iRettig fc Pohll 
I2012t IRessler et al.ll2014) , and this mechanism is the only 
possible one if the rims are equally thin at radio wave¬ 
lengths, as is true in some regions of Tycho. Thus, the 
age-limited interpretation of our NuSTAR results would 
demand low magnetic field strengths and magnetic-field 
damping to be responsible for the X-ray rim morphology. 
Howe ver, a recent detailed analysis of Tycho (jTran et al.l 
l2015fl finds that even in damping models, magnetic field 
strengths /^30 fiG are required to fit the radio and X-ray 
profiles of the thin rims. 

Alternatively, the electrons are loss limited, and the 
apparent trend of Aroiioff with Vs arises because of obliq¬ 
uity effects. In the scalings relating Aroiioff to Vg, we have 
assumed constant Rj, yet the shock obliquity 0Bn is ex¬ 


pected to change around the periphery of the SNR as the 
shock encounters a roughl y uniform magn etic field in the 
environs of a Type la SN (jRevnoldill998f) . For example, 
changes in obliquity are thought to be the origin of the 
X-ray bright, synchrotron rims aligne d in the northeast 
and southwest directions o f SN 1006 (jRothenfing et al.l 
12004 iRevnoso et al.l 1201311 . The magnetic field is not 
as ordered in Tych o as in SN 1006 (see Figure 8 of 
IRevnoso et al.lll9T^ . and further work is necessary to 
assess the obliquity dependence of the particle accelera¬ 
tion in Tycho. 


5. CONCLUSIONS 

We have reported the large, ^750 ks NuSTAR observ¬ 
ing program toward the historical supernova remnant 
Tycho. Using these data, we produced narrow-band im¬ 
ages over several energy bands to identify the locations 
of the hardest X-rays and to search for radioactive de¬ 
cay line emission from '^"‘Ti. We find that the hardest 
>10 keV X-rays are concentrated to the southwest of 
Tycho, where recent Chandra observations have revealed 
high emissivity “stripes” associated with particles ac¬ 
celerated to the knee of the cosmic-ray spectrum. We 
do not find evidence of “^^Ti, and we set limits on the 
presence of “^^Ti, depending on the velocity and distri¬ 
bution of the metal. In order to be consistent with the 
reported Swift/BAT detection, the Ti must be expanding 
at moderate-to-high velocities and/or distributed over 
the majority of the SNR. Our spatially-resolved spectro¬ 
scopic analyses of sixty-six regions showed that the high¬ 
est energy electrons are accelerated by the fastest shocks 
and in the lowest density regions of the SNR. We find 
a steep dependence of the roll-off frequency with shock 
velocity that is consistent with the maximum energy of 
accelerated electrons being limited by the age of the SNR 
rather than by synchrotron losses, contrary to previous 
results obtained for Tycho. One way to reconcile these 
discrepant findings is through shock obliquity effects, and 
future observational work is necessary to explore the role 
of obliquity in the particle acceleration process. 
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